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1 Introduction 

Since the conditions for generating a matter-antimatter asymmetry were laid out by Sakharov 
[1], many scenarios for baryogenesis have been proposed, based on high-energy processes in 
the early Universe (for a review, see for instance [2]). Here we report on progress of a tech- 
nique that can be used to examine various models of electroweak-scale baryogenesis [3-5], 
whose common feature is that the source of baryon number violation is the electroweak 
anomaly [6]. 

All the key symmetry violations (CP, C and baryon number) occur in the Standard 
Model, and the requirement of departure from thermal equilibrium is determined by early- 
Universe dynamics and the electroweak symmetry breaking transition. In the Standard 
Model (SM), CP violation is provided by the CKM matrix of quark mixings [7, 8] and 
possibly by the PMNS matrix of neutrino mixings [9, 10], although this is less well mea- 
sured experimentally. Beyond the SM, CP-violation can be provided for instance through 
Major ana masses, multiple Higgs field with complex coupling and/or higher-dimensional 
effective terms. 

In electroweak baryogenesis the baryon asymmetry is expected to occur at the elec- 
troweak scale of ~100GeV, at which time the Hubble expansion rate is around 10~^eV 
- very slow by the standards of microscopic electroweak processes. Given that a ther- 
mal electroweak phase transition in the Standard Model cannot be first order for current 
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Higgs-mass bounds [11, 12] , we are left with the problem of finding dynamics that can 
cause sufficient departure from thermal equilibrium. Such a possibility is raised by models 
inspired by preheating at the end of inflation, where the phase transition is non-thermal 
[13-20], or in extensions of the Standard Model with a more complicated scalar sector, 
where the phase transition can be strongly first order [21]. 

Whilst it is nowadays a relatively straightforward task to simulate the classical dy- 
namics of the bosonic degrees of freedom of the Standard Model and minimal extensions, 
using techniques from lattice gauge theory, the fermions are more problematic. And given 
that all the CP violation of the Standard Model is in the fermion sector, and the baryon 
number itself is a fermion quantum number these cannot merely be ignored. One approach 
is to "integrate out" the fermions, and treat their dynamics as new terms in an effective 
action. This approach has been taken in a number of studies [22-25] (see also [26-28]). 

Here we present a numerical method that evolves the fermions along with the bosons 
and, when extended to three families of fermions, can include CP violation from the fermion 
sector directly. 

The technique that we use to simulate the dynamical fermions was first laid out in 
[29], with a sample study of the effect of fermions on the evolution of oscillons in 2+1 
dimensions. The idea is to model the quantum averages ( |...| ) by ensemble averages (...)e, 
and so by simulating an ensemble of fermion realizations, one is effectively sampling the 
evolution of the whole set of fermion mode functions. The utility of this is that the number 
of mode functions varies as in three-dimensions (here Ux is the number of lattice sites in 
one direction), which is prohibitively expensive numerically, whereas the number of modes 
required to give a reasonable sample is expected to vary as [29] (see also the 3+ ID 
simulations in a scalar-fermion theory in [30]). 

One check that the method is viable for reproducing the quantum anomaly has al- 
ready been performed in an axial Yang-Mills-Higgs-fermion system in 1+1 dimensions [31], 
where the baryon-number violating processes where directly simulated. The results of the 
ensemble method were in complete agreement with the evolution using the full set of mode 
functions [32, 33]. While there is no numerical gain in 1+1 dimensions by using the ensem- 
ble method, the results of [31] gives confidence in the technique, and provided motivation 
for the present 3+1 study, which would not be practical using the full set of mode functions, 
certainly for lattices of the size used in [24]. 

In this paper we use the fermion ensemble method to simulate the electroweak quantum 
anomaly directly. In this first study we shall be using only a single family of Standard Model 
fermions ^ which is enough to establish the viability of the method by direct reproduction 
of the anomaly. We shall also be omitting the U(l) hypercharge and SU(3) colour of the 
Standard Model, as this does not affect the anomaly process. 

The structure of the paper is as follows: We start in section 2 with a presentation of 
the continuum model that we aim to simulate and in section 3 discuss how we treat the 
fermions and bosons numerically. In an equilibrium environment, the key baryon number 
violating process is the sphaleron transition. In section 4 we test the real-time fermions in 

^We include a right-handed neutrino in our definition of the Standard Model. 
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the background of handmade sphalerons. We then describe how to do full boson-fermion 
dynamics in the context of a fast electroweak transition in 5, where we demonstrate the 
approach with sample simulations. We conclude in section 6. In order to keep a flow to 
the paper we have included a number of the technical details in the appendices, including 
a list of conventions in App. A, the lattice implementation in App. B, the fermion initial 
conditons in App. C and the details of the handmade sphalerons in App. D 

2 The SU(2)-Higgs model with chiral fermions in 3+1 dimensions 

Our field content is one family of the Standard Model (extended to include a right-handed 
neutrino), but without U(l) hypercharge, SU(3) gluons and without colour degrees of 
freedom. Therefore we have a complex Higgs doublet 0; SU(2) gauge fields Wj^; a left- 
handed SU(2)-doublet quark field qi = (u^jdi); two right-handed SU(2)-singlet quark 
fields ur, dfi; a left-handed SU(2)-doublet lepton field II = {i^Li^l}', and two right-handed 
SU(2)-singlet lepton fields e/j, vr. The continuum action is then written as 

S = Sh + Sw + Sf + Sy, (2.1) 
where the set of different components is given by 



Sh = — I d X 



(2.2) 



Sw = - I d^x -WH^W'^"", (2.3) 



1, 
4 

Sf = - j d'^x [qLl^D^qL + UR^i^D^UR + dR-fi'Di.dR 

+Ili''D^Il + VRl^D^UR + eR^^'D^CR] , (2.4) 

Sy = - J d^x [c^'qLct'UR + G^'qict'dR + G%ct)eR + G%ct)UR (2.5) 

+G''qL^UR + G'^qL^R + G%^eR + GnL4>i^R 
+h.c.] . 

The absence of a U(l) hypercharge means that one is allowed more Yukawa terms than 
the equivalent Standard Model action, but we simply set to zero those that would not be 
allowed had we included the Standard Model U(l). 

The charges of the fields are determined by the covariant derivatives, which are 

D,<P = (d, - fa'^W^;) (2.6) 

Df^qL = (df,-'-^a^W^^qL, D^ur = O^^ur, D^dR = d^dR, (2.7) 

D^^Il = (d^ - y^r'VF;^ II D^eR = d^eR, D^vr = O^ur, (2.8) 
and the SU(2) field-strength is defined by 

[D^,D,]<^ = -|a'^VF>. (2.9) 
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As well as the conserved current following from the presence of the gauge symmetry, there is 
a baryon current and a lepton current coming from the global symmetries q — ?• exp(ia7^)(; 
and / — )• exp(iQ;7^)/, 



3(h) 

Hi) 



(2.10) 
(2.11) 



Classically, these currents are conserved. Quantum mechanically however, one finds that 
[34] 



d, 



32^ 



-f w°- w°- 



where 



167r2 



-pupa 



(2.12) 
(2.13) 

(2.14) 



and is the number of fermion families, which we have equal to one for our simulation. 
We then find that the baryon number, Nf = J d^xj^^^y is related to the Chern-Simons 
number, Ncs = / d'^xK^, as follows 



N 



f 



Ncs- 



(2.15) 



The main aim of this paper is to confirm that this truly quantum mechanical relation holds 
in our numerical simulations, and to apply the method to a fast electroweak transition. 
We will also calculate the average Higgs field 



1 

Vol 



d^x(h\ 



and the Higgs winding number, 

Nw = ^ J d'x e,jkV^d,VV^d,VV^dkV, V = ^ 



(2.16) 



(2.17) 



For convenience: Partial charge conjugation transformation 

On a lattice it is is more convenient to redefine the Fermi fields so that their kinetic terms 
are just the standard Dirac form, i.e. the fermion-gauge interactions are vector-like rather 
than chiral. This may be achieved by the following definitions [32, 33, 35, 36], 



*L 

XR 
XL 
iR 



UR, 

dR, 



II 



r>-l -T 
eR = C Xl 



^R 



dR 



^R 



-xlc, 



(2.18) 
(2.19) 
(2.20) 
(2.21) 
(2.22) 
(2.23) 
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We then find that 



(2.24) 
(2.25) 
(2.26) 



iLl^a'^W^h = ^Ri^a'^W^-^R, (2.27) 
where some integration by parts has been done on the first three equations, leaving us with 



Sf 
Sy 



(2.28) 



C^^c^PrC + G^x(i>^PR^^ + C'^^Prx - G''i4)^PR^ + h.c] . (2.29) 



The end result is that whereas we before had two left-handed doublets and four right- 
handed singlets, these are now collected into one full Dirac doublet and two Dirac singlets. 
The singlets only interact via the Yukawa term. It is straightforward to check that in this 
formulation one has 



'(5) 



C— conjugated 



Original 



(2.30) 



3 Modelling the bosons and fermions 

The continuum model above is discretized on the lattice as described in App. B. For 
the bosonic fields the classical Hamiltonian equations are derived, with the understanding 
that fermion bilinears are to be represented by quantum averages. The fermion evolution 
equations are also found in the usual way by variation of the action, and amount to a linear 
Dirac equation in a time- and space-dependent bosonic background. 

One issue to address are the lattice fermion doublers. For every physical fermion mode, 
the standard lattice discretization prescription generates a set of fifteen doublers, which 
count as real degrees of freedom at finite lattice spacing [37] ; in the context of the anomaly 
these are particularly troublesome. Since they contribute with opposite signs, the doublers 
cancel the anomaly from the physical fermions, and set the total anomaly to zero [38] (see 
also Fig. 3). 

Fortunately, a simple approach exists to breaking the symmetry between these dou- 
blers, namely adding a Wilson term to the lattice action (B.13). As was the case in 1+1 
dimensions [31-33], we find that it is sufficient to add a spatial Wilson term to cancel the 
space-like doublers, and to not initialize the time-like doublers. For small enough timestep 
in the numerical evolution, the time-like doublers stay unexcited, at least for the duration 
of our simulations (see also section 4). 

The next step is to replace quantum averages by ensemble averages, as discussed 
in detail in [29, 31]. The process starts by noting that the continuum canonical anti- 
commutation relations for continuum fermions are 

{'0a(t,^),V'fl(i,^')} = ^a/3'5(^-^')i (3-1) 



- 5 - 



and we may expand the field operator as 

dPp 1 



■ip.x 



in terms of tlie creation and annihilation operators. The anti-commutation relation are 
then equivalent to imposing 

{hl{p)M^)} = {27rf{2ujp)6rsS{p-p;), (3.2) 
= {27rf{2ujp)SrsS{p-p;), (3.3) 

and we may calculate bi-linears such as 

where fi is the fermion mass; it is these kinds of quantities that we aim to reproduce using 
ensemble averages. 

To model the quantum averages we introduce two ensembles of fermions, M(ale) and 
F(emale), according to the mode expansion 

VM,F(t, x) = ^T.j ^ [Up)UsW^-^ ± r/,(^)K(p)e-^-] , (3.6) 

and where the exact same random numbers ^, C, are used in a given Male and Female pair. 
By then requiring that the variables ^ and rj satisfy the ensemble average relations 

{Up)es{^))e = {2T,f{2u:p)5rs5{p - p'), (3.7) 
{rir{p)r]W))e = {27rf{2iOp)6rsS{p-p;), (3.8) 

(where 6{p — p') is understood to refer to the lattice version of the Dirac delta), we may 
for example calculate 

(2^)3 2w/ ' ^^-^^ 

which leads to the simple prescription 

( mxmy)\ ) ^ _ (^M(x)^Hy) + (MoF)). _ 

For more details of the method see [29] . 

The mode expansion described above constitutes the initial conditions for our fermions, 
which then amounts to generating sets of and inserting them into eq. (3.6). In the 
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Figure 1: The fermion number compared to the Chern-Simons number in a handmade 
sphaleron transition, for different sizes of the fermion ensemble, Nq. There is clear conver- 
gence to the exact result, which is within statistical error bars (not shown) for all values 
of the ensemble here, amn = 0.42. 



case of a general Yukawa interaction and CKM mixing matrix, one has to take care to 
diagonalize the fermions first to initialize the mass eigenmodes; this is explained in more 
detail in App. C. For the simulations presented here, we will restrict ourselves to Yukawa 
couplings proportional to the identity (see also App. B) 

G*^ = = G"^ = -G" = Ayuk- (3.11) 
4 Sphaleron transitions 

Now that we have the full prescription for evolving the system of bosons and fermions, we 
wish to test whether the method is reliable for the anomaly dynamics driven by sphaleron 
processes. We are mainly interested in whether the numerical fermions react in the correct 
way as the SU(2) gauge field undergoes a change in winding. To that end we shall start by 
evolving the Higgs and gauge fields by hand, rather than through their equations of motion, 
so we need to know how to construct a winding event. This is the method used in 1+1 
dimensions [31-33], and here we extend the construction to the SU(2) sphaleron following 
the techniques of [39] (see also [40] for a nice exposition). The detailed implementation 
can be found in App. D. 

The end result is a family of configurations starting and ending in a pure-gauge vacuum 
of the gauge-Higgs theory, which takes the fields through a sphaleron transition. In this 
way there is a relative winding number between the vacua, both in terms of winding in the 
Higgs field, and the Chern-Simons number of the gauge field. On this bosonic background 
we evolve the fermions according to their dynamical evolution equations, and test that 
they respond in the expected way and acquire a non-zero fermion number in accordance 
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Figure 2: Handmade sphaleron transitions for different values of the lattice spacing 
aniH = 0.105 — 0.42. The lattice spacing dependence is well under control. The lattice size 
was kept constant at rix = 32. The band is the statistical errors. 



with the anomaly equation. The fermions are initialized in their vacuum as described in 
App. C. 

For the moment, we will set the Yukawa couplings to zero. We generate a double- 
sphaleron trajectory on a = 32'^ lattice with niH/mw = 2, a timestep at/ux = 0.05 
a lattice spacing amn = 0.42 and a Wilson coefficient r = 1.0. Fig. 1 shows Chern- 
Simons number A'cs in time (black) and the combination (B + L)/2 increasing the fermion 
ensemble. The coloured lines show the ensemble averaged anomaly for Nq = 640, 1280, 2560 
and our largest ensemble Nq = 10240. We see that the anomaly is very well reproduced 
at the ten percent level for these ensemble sizes, and that this can be improved upon 
in a straightforward way. For a 32'^ lattice, numerical efficiency is only gained over the 
"non-statistical" approach for ensemble sizes less that 32^ = 32768. 

The next test is to vary the lattice spacing (for fixed number of lattice points) , shown 
in Fig. 2. We see that the anomaly is well reproduced for all lattice spacings, although the 
am-H = 0.105 is a little low, possibly because the physical volume is quite small UxamH = 
3.36. Also the largest lattice spacing amn = 0.42 is about 5 percent low, which may in 
turn be put down to discretization errors. We note that the anomaly is expected to be one 
of the "hardest" observables to get right, as it does not involve any volume averaging (such 
as for instance particle number in a homogeneous background would [30]). Also, as can be 
deduced from the level-crossing picture of the anomaly, all momentum modes contribute to 
the anomaly, both in the IR and the UV. Hence the statistical averaging has to be accurate 
for all modes independently, to get the correct anomaly out. It is not possible to leave out 
some of the UV modes, as if one were only interested in the IR physics. In this light, we 
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Figure 3: The anomaly as the spatial Wilson coefficient is varied. 

consider it remarkable that our simulations reproduce the anomaly so well. 

In Fig. 3 we show how the anomaly depends on the coefficient of the Wilson term. As 
mentioned, for r = lattice doublers are expected to pair up and cancel out the anomaly, 
whereas for r = 1, the anomaly is manifest. We indeed see that the anomaly is identically 
zero for r = 0, but also that once r > 0.25 the anomaly is largely unaffected by r (within 
statistical errors). This is as we would have expected. Also note that the time-like doublers 
are not cancelled by a Wilson term, but that these are not excited enough to cause problems. 

We then consider the case where scalars are directly coupled to the fermions, i.e. the 
fermions are massive. As shown in Fig. 4, at finite Yukawa coupling convergence slows down 
with increasing coupling. This is completely analogous to the situation in 1+1 dimensions 
[31], and the solution is also the same: increase the size of the ensemble. Nq = 10240 is 
again sufficient for Ayuk < 0.1, corresponding to a fermion mass of Ayukt'/\/2 ~ 17GeV, 
i.e. all Standard Model quarks except the top. But to include the top, Ayuk — !> we 
need many more realizations. We did a test at the smaller lattice size of = 16^ with 
Nq = 8 X 10240 and indeed found that the anomaly is still present, shown in Fig. 5. Such 
simulations suffer from having small physical volume, and the fermion number is a little 
low. We also explicitly checked that decreasing the timestep and/or increasing the Wilson 
coefficient r did not improve convergence. [39] We can conclude from these tests that the 
fermion ensemble method works even for very subtle quantum observables like the anomaly, 
and that we need Nq = 0(10"^) to get a measurement at the ten percent accuracy level, at 
zero and small Yukawa couplings, up to a fermion mass of a few GeV. A larger ensemble 
is required to consistently include heavy quarks, i.e. the top. Lattice issues like lattice 
spacing dependence and fermion doublers are well under control. We will now proceed to 
also make the bosonic fields dynamical and include the fermion back-reaction. 
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Figure 4: Handmade sphaleron simulations, varying the Yukawa coupling. Nq = 10240. 
Ux = 32, aniH = 0.42. The Yukawa couplings correspond to fermion masses of m/ = 
87, 44, 22, 11 GeV, respectively. 



5 Including fermion backreaction: tachyonic electroweak transition 

In the tests with handmade sphalerons presented above, the fermions were just "along for 
the ride", to see how well the fermion equations of motion behave in a time-dependent 
bosonic background. But for baryogenesis, we need the bosonic fields to be dynamical, 
too, and for the fermions to back-react. Indeed, in the SM CP-violation is a fermion 
back-reaction effect. 

As a testing ground for this, we consider a cold tachyonic electroweak transition, where 
the Higgs field experiences a fast quench. The Higgs potential flips instantaneously as 

V{cP)=+f,^U, (t<0) ^ AUV-!i-j , (t>0). (5.1) 

The Higgs field is generated by Monte-Carlo sampling of a classical ensemble, reproducing 
the quantum correlation functions at zero temperature in the symmetric phase [18, 19, 23]. 
The fermion fields are generated as an ensemble of Nq vacuum realizations at = 0, as 
above. The initial gauge fields are then found by setting A^{t = 0) = 0, and solving the 
lattice Gauss law in the background of the generated Higgs and averaged fermion fields. 

The dynamics start at the moment when the Higgs potential is flipped, and the low- 
momentum modes k < ^ experience a tachyonic (or spinodal) instability, and start growing 
exponentially. This leads to a strongly out-of equilibrium electroweak transition, where 
large particle numbers, effective diffusion of Chern-Simons number and in the presence of 
CP- violation, baryogenesis can occur [13, 20, 23, 24]. In the present work, we have only 
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implemented one generation of Standard Model fermions, and so we are not able to include 
CP-violation via the CKM matrix. We will leave that for future work. 

As in [23] we will monitor the volume-averaged Higgs field, the Chern-Simons number 
and the Higgs winding number. But rather than infer the baryon and lepton number 
through the anomaly equation, we now calculate it directly. We note that the anomaly 
equation states that changes in the Chern-Simons and Baryon/Lepton number obey 

AB = AL = AN,, ^ AiVes = ^^ + ^^ _ (5.2) 

Fig. 6 shows a tachyonic run where the fermions do not back-react (dashed) and one 
where they do (full). Yukawa couplings are zero, and Nq = 5120. We see that the Higgs 
field (black) "falls off the hill" and oscillates around its broken phase minimum. Meanwhile, 
the Chern-Simons number (red) and Higgs winding number (green) bounce around out of 
equilibrium, but eventually settle down in a common minimum. To check our numerics, 
we also monitored Gauss law, which throughout is conserved to computer accuracy. As for 
the sphalerons above, we see that the fermions (blue) reproduce the anomaly rather well, 
however, the gauge and Higgs fields are now dynamical. 

What is particularly interesting is the effect of having fermions affecting the bosonic 
dynamics. At early times, the evolution of the bosonic observables is unchanged. But then 
around rriHt = 15, the trajectories of both fermion number and Chern-Simons number di- 
verge, depending on whether the bosons feel the fermions, eventually ending up at different 
minima, and hence a different number of fermions. 

Finally, in Fig. 7 we show a run with small enough Ayuk = 1/128, that Nq = 10240 gives 
reliable convergence. This corresponds to a 1.4 GeV fermion, i.e. slightly heavier than the 
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Figure 6: A tachyonic transition with (full) and without (dashed) fermion back-reaction, 
starting from the exact same initial condition. Shown are the trajectories of the Chern- 
Simons number (red), the average Higgs field squared (black) the Higgs winding number 
(green) and the fermion number (blue). 




Figure 7: A full tachyonic run including fermion backreaction and Yukawa coupling. The 
fermion mass is m/ ~ 1.4 GeV. 

charm quark. As for the sphaleron runs, we see that the agreement is only approximate, 
but clearly the anomaly is reproduced, improveable with larger statistics. 
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6 Conclusions 



While there have been many simulations of the bosonic fields in the Standard Model, this 
is the first that has included the fermionic sector, rather than integrating them out and/or 
treating them as passive observers to the gauge-scalar dynamics. This has been achieved 
by applying the technique pioneered by Borsanyi and Hindmarsh [29], where the many 
mode-functions of the fermion fields, required for a full quantum treatment, are replaced 
by an average over a statistical ensemble of Male and Female fields. 

We have performed a detailed study of the lattice parameters, and found that the 
fermion number does indeed reproduce the quantum anomaly, with the fermion number 
tracking the Chern-Simons number through sphaleron transitions. We also simulated a 
fully dynamical fast electroweak transition including fermions, and found that the back- 
reaction of fermions can indeed significantly alter the evolution of the gauge fields. 

From a technical point of view, the numerical effort is significant, in that implementing 
Nq fermion realizations on a n^. lattice requires x Ng x 1.6 kB, so that for our single 
generation, no-colour full run of Fig. 6 we use approximately 270 GB. It is however impor- 
tant to note that the anomaly is probably one of the most sensitive observables, since it 
requires every single lattice mode to be very precisely reproduced. This means that for 
studying the (IR) dynamics of electroweak baryogenesis, we can probably make do with 
a somewhat smaller Nq. Alternatively, one may consider returning to the non-ensemble 
approach [32, 33], where one is however hampered by the need for large lattices to correctly 
simulating tachyonic transitions. 

The next step is now to study dynamics both for baryogenesis and electroweak preheat- 
ing; including all three generations, mixing and CP-violation for both the quark (CKM) 
and lepton (PMNS) sector as relevant for the most minimal version of Cold Electroweak 
Baryogenesis; but also for electroweak baryogenesis in more general setups like Standard 
Model + second Higgs doublet or Standard Model -|- scalar singlet. Demonstrating that 
this is a realistic aim (also in practice) is perhaps the most significant result of the present 
work. 
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A Conventions 

We consider 3-l-lD Minkowski space, with the metric (— , -|-, -|-, -|-), and Dirac algebra 

{7^7'^} = 2r?^^ ^ = iVV, 75 = -i7V7V, (A.l) 
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with in the Weyl representation 



C=(~f .\]= -nV, CCt = 1, Ct = -C, (A.3) 
\ la I 

e = ^^2, e-Ve = -^r"*, (A.5) 

V'R = ^(1 + 7')V', V'L = ^(1-7')V'. (A.6) 



The charge-conjugate scalar doublet is given by the complex conjugate 



B Lattice action 



(A.7) 



The theory is discretized on a lattice with spatial lattice spacing ax, timelike spacing 
at = dtttx- We find it convenient to work with the following rescaled and re-defined variables 

^lattice = ^W^continuum , 

*<^' - ( ) - ^ ^« ( % ) - • 

^(x) ^ a^/^^{x), x{x) ^ a3/2x(x), ^ a3/2e(x), (B.3) 

thus making them dimensionless. Then we have the link variables, 

Uo{x) = exp{-iata''WS{x)/2), (B.4) 
Ui{x) = exp(-m^a"W,"(x)/2), (B.5) 

Uij{x) = Ui{x)Uj{x + i)uj{x + j)u]{x) = exp {-i]^a^ a^W^^ + ...j (B.6) 

and corresponding electric fields 

Ei{x) = -Ui{x)ul{x + Qi), Ef{x) = Tria''Ei{x), (B.7) 
= -Ujix + 0)C/,(a;), (x) = Tr ia^Eiix). (B.8) 

We need the derivatives, where we will not include the lattice spacings (they are compiled 
into lattice coefficients, see below) 

7r{x) = <^{x + 0)-<^{x), (B.9) 
D^^ = U^{x)^{x + /u) - ^{x), 
D'^<^ = $(x)-C/t(x-/i)$(x-/x), 
D'^Df,^ = U^{x)^{x + 2^>(x) + U^x - n)^x - fi), 

D, = ^[D, + D'^]. (B.IO) 
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Note that we are using the compact formulation of the gauge action, and have chosen to 
evolve using temporal gauge Wq = {Uq = I). In order to remove the spatial fermion 
doubler modes we will employ Wilson fermions in space, for which we need the Wilson 
term [32, 33, 41] 

^W^ = -^^^D',Di-i', ^Wx = -^^xd'Ax, ^WC = -l^^d'M, (B.ll) 
where r is a constant parameter. We will also introduce the quantities 



at a 



la lat {amnf ^H^m 

T — > Ph = PR = T — : -A = ^ {ti.l2) 



The action is then given by 

S = Sh + Sa + Sp + Sw + Sy, 

with 



TV 



x,t 



(B.13) 



2\ 2 



(B.14) 
(B.15) 



Sf = Y.- [^7°^o^ + X7°^oX + ^Y^oe] - ^ [^t'A^ + Xl'~diX + h'd.i] , (B.16) 



x,t 



(B.17) 



To which we add the Yukawa terms, now written in terms of the matrix <I>, and where we 
for simplicity impose that 



and in addition set 
Then we have, using 

the addition to the action 



yuk' 



/3y 



at 1 
a ^/x' 



(B.18) 



(B.19) 



(B.20) 



Syuk = - E 0^ + ix, Ob^lj^) ■ (B.21 



The equations of motion follow by variation. For the Higgs field 



-Tr$T<l) 



Ph 



(B.22) 
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where 



and 

are ensemble average over fermion realizations. For the gauge field we have 



(B.23) 

(B.24) 
(B.25) 



+ 



rat 



„ L J apg L 



a/3* L 



^yia'^Uy.n'if 



y,n ^ y+n 



^y+nUlja'^^y 



(B.26) 



and we have Gauss law 



2^ 

Ph Pg 



Trier" 



[^'y^^'ia'^^'y+Q + ^y+Q-i'^ia'^^'y] . (B.27) 



Finally, the fermion fields evolve according to 



1 



[dt + d',]^a 



1 



[dt + d[]ix,0'' 



2^-7°^ [A + D'^^a + ^7°A'A^a - (3y^abl\x,0l (B.28) 



2' " 2a' ' ' ' 2a 

The currents that are particularly relevant to us are the baryon and lepton currents 
which, when converted to male and female fermions, becomes 



4) 



-- [-^M(x)7°7'^'F(a; + 0) - ^m{x + 0)7°7'^F(a;: 

+Xm(2;)7°7^Xf(x + 0) + xm{x + 0)7°7^xf(x) 
+^Mixh^l'CF{x + 0) + ^m{x + 0)7°7'eF(x) 
-f'F(a;)7V^Af(a^ + 0) - ^f{x + 0)7V^a/(2;) 
+Xf(2;)7Vxm(x + 0) + xf{x + 0)^^-f^XM{x) 

+e>(x)7%5eM(x + 0) + + o)7%'eM(x)] , 



(B.30) 



J(5) 



1^ 

'4 a 



-^'M(2;)7''7'f^fc(a;)^F(x + k)- ^m{x + k)-f''-f^ul{x)^F{x) (B.31) 

+Xm(x)7N^Xf(x + A;) + XAf(x + /c)7S^Xf(2;) 
+^M{xh''l^^F{x + k)+ ^m{x + k)-f''-f^^F{x) 
-^F(x)7S^f/fe(x)^'M(x + k)- ^f{x + k)j''-r^ul{x)^M{x) 
+XF{xh''j^XM{x + k) + XF{x + k)j''j^XM{x) 
+^f{x)j''j^U.i{x + k) + ^f{x + k)j''-r^^M{x) 
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C Fermion initial conditions 



In order to specify the initial conditions we need some mode functions for the fermion 
fields. The ones we use are those that are vacuum fluctuations for a Higgs field taking the 
value 4) = (0, f/v^), 4> = ("y/y^, 0). The lattice equations of motion then read 



X \ vat 



V2 U r 



(C.l) 

(C.2) 
(C.3) 



with M and M some general Yukawa mixing matrices. This is a set of linear equations, so 
we perform a Fourier decomposition and look at positive frequency modes 

^+ = U{k)e'^-'', (C.4) 

and find that by defining 

So = sin(A;o at) /at = - sm.{uj at) /at, Si{k) = sm{ki)/a, m{k) = r ^(1 - cos(/ci)), 

(C.5) 



one gets 



a V2 \ U^^> 

' \ _ m{k) ( U^^^ 



(C.6) 



(C.7) 



now define fh = m{k)/a and v = v/\pl and make the simplification that the Yukawa 
couplings are related by 



-G° 



so that defining 



(C.8) 



(C.9) 



allows us to diagonalize the equations by introducing the canonically normalized Fermi 
fields Oj, 

1 



9.1 



1 
1 
1 

7%g^ 



(C.IO) 
(C.ll) 
(C.12) 
(C.13) 
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The positive frequency modes of these Fermi fields then satisfy 

-irs^.Ul = {m + Wf)U^, (C.14) 

-i7%f/2 = (m- |5l2)C/2, (C.15) 

-i7%f/3 = (?^+bT)[/3, (C.16) 

-i7%f/4 = (m - |5T)C/4. (C.17) 

They are now in the usual form, and we can solve for the canonical mode functions, where 
in general 

{-i-i^'s^- ^l)U{k) = Q, (C.18) 

{i^s^- li)V{k) = (C.19) 

fiO = ^./IflT^^ (C.20) 



which are 



and i]s 




D Handmade Sphalerons 

The basic idea [39] is that the vacuum manifold of the Higgs scalar is a three-sphere, S^ac) 
and one wants to construct some other three-sphere in order to have non-trivial maps 
between them. At spatial infinity the Higgs field is in the vacuum, but spatial infinity is 
only a two-sphere, S^, and the maps S'^ — y ^vac '^^^ topologically trivial. However we can 
construct a three-sphere by considering a loop in configuration space, Sf^^, and forming 
the smash product 5^ A Sf^^^, which is a three-sphere, and so one may construct non- 
trivial maps. This loop in configuration space means we have a family of configurations 
parametrized by an angle. 

In practise this may done by parametrizing the S'^ A Sf^^ three-sphere using the unit 
four-vector 

Pa = (sin r sin 9 cos ip, sin T sin sin (p, cos^ T + sin^ F cos 9, sin F cos F(cos 9 — 1)), 

which is an explicit representation of the smash product, where F parametrizes Sf^^^, and 
9, ip parametrize S"^. ^ 

Now we work with the usual complex doublet and write the set of configurations at 
infinity as 



V I sin F sin 9e^^ 

\ e~^^ (cos F + i sin F cos ( 




(D.l) 



^Note that 0<r<7r, 0<i^< 2-k, < 6 < tv and that the circle at S = is mapped to the point 
(0,0,1,0), and the two-sphere at F = is mapped to the same point (0,0,1,0), thus confirming the smash 
product nature of the parametrization. 
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which leads to the ansatz for the fuh radial dependence of the family of configurations 

Mr, 0, = /i(r)<Aoo(r, e,^) + ^[l- h{r)] ^_ J^^^ ^ ^ , 

where the radial profile function is a smooth monotonic function obeying h{0) = 0, h{oo) = 
1. This leads to 

Mo,n{r,e,ip) = ^i^] , (D.2) 



V2\ 



sm 



cl)^=^/2{r,e,^) = -^h{r){ ), (D.3) 



cos 



so we see that the start and end points of the F loop just correspond to the vacuum, with 
the centre of the loop (T = 7r/2) being the sphaleron [39]. 

Of course, we also need an ansatz for the gauge field, and this is achieved by construct- 
ing the unitary matrix 

n^iT,e,^) = ^( ^l- (D.4) 

The gauge field is in the vacuum at infinity, and so we take the asymptotic ansatz to be 

TF,(oo)(r,0,(^) = -'-dinooiT,9,^)nl{T,e,^), (d.5) 

with the full radial dependence of the family of gauge fields being given by 

Wj{r,r,e,ip) = /(r)VF,(^)(r,r,0,<^), (D.6) 
where /(r) is a smooth monotonic function satisfying /(O) = 0, /(oo) = 1. 
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